Bulk Viscosity in Hybrid Stars 
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We compute the bulk viscosity of a mixed quark-hadron phase. In the first scenario to be discussed, 
the mixed phase occurs at large densities and we assume that it is composed of a mixing of hyperonic 
matter and quarks in the Color Flavor Locked phase. In a second scenario, the mixed phase occurs 
at lower densities and it is composed of a mixing of nucleons and unpaired quark matter. We have 
also investigated the effect of a non- vanishing surface tension at the interface between hadronic and 
quark matter. In both scenarios, the bulk viscosity is large when the surface tension is absent, while 
the value of the viscosity reduces in the second scenario when a finite value for the surface tension 
is taken into account. In all cases, the r-mode instabilities of the corresponding hybrid star are 
suppressed. 

PACS numbers: 97.60.Jd, 26.60.+C, 25.75.Nq, 04.30.Dg 



I. INTRODUCTION 

The discovery by Andersson, Friedman and Morsink 
of r-mode instabilities in neutron stars put rather severe 
limits on the highest rotation frequency of pulsars Q, • 
These constraints can be incompatible with the existence 
of millisecond pulsars, if the instability is not suppressed 
by a sufficiently large viscosity. Actually, for a star com- 
posed only of neutrons and protons and for temperatures 
larger than roughly 10 10 K, the bulk viscosity due to the 
modified Urea process is large enough to damp the insta- 
bility 0. On the other hand, when the star cools down 
to lower temperatures, the instability is not suppressed 
and the star is forced to loose angular momentum via 
emission of gravitational waves Q . More recently, it has 
been noticed that the existence of viscous boundary lay- 
ers at the interface between the fluid core and the crust 
can stabilize the star [2jj if the rotation period is longer 
than ~ 1.5 ms 0, At very low temperatures, below 
10 8 K, shear viscosity becomes large and it allows older 
stars to increase their angular velocity by mass accretion. 

Actually, compact stars can be constituted by a larger 
variety of particles than just neutrons and protons. One 
possibility is that hyperons form at the center of the star. 
It has been shown that, due to non-leptonic weak reac- 
tionSjbulk viscosity can be rather large for an hyperonic 
star 0, @ , which therefore can emit gravitational waves 
only if its temperature is ~ 10 10 K and its frequency is 
larger than 10-30% of its Keplerian frequency. 

The formation of quark matter inside a compact star 
has been discussed extensively in the literature. Bulk 
viscosity of non-interacting strange quark matter is very 
large On the other hand, recent studies taking 

into account quark-quark interaction revealed the possi- 
ble existence of color superconducting phases in which 
quarks form Cooper pairs with gaps as large as 100 
MeV 12]. In that case, bulk viscosity is strongly sup- 
pressed by the large energy gaps, r-mode instabilities are 
not damped and pure color-superconducting quark stars 
seem therefore to be ruled out by the pulsar data (l3|. 

In the present paper we are interested in studying the 



viscosity and stability, respect to r- modes, of an Hybrid 
Star (HyS), for which no quantitative analysis has been 
performed so far. If the star is made of hyperonic matter 
and of non- interacting quarks, it is rather obvious that 
the viscosity of the HyS should be large, due to the large 
viscosity of its constituents. The real question concerns 
the case in which color-superconducting quark matter is 
considered, since its viscosity is negligible. There are in 
principle (at least) two sources of viscosity for HySs. One 
originates at the interface between the crust and the fluid 
interior [13j and it is similar to the one discussed above 
in the case of purely hadronic stars 0,0. In our work we 
will not discuss that possibility and we will instead con- 
centrate on a quantitative evaluation of the bulk viscosity 
in a mixed phase (MP) composed either of hyperons and 
color-superconducting quark matter or of nucleons and 
non-interacting quarks. We have also investigated the 
effect of the existence of finite-size structures in the MP. 
These geometrical structures can exist if the surface ten- 
sion at the interface between hadronic and quark matter 
is not vanishing. We will show that the bulk viscosity of 
the MP is in general rather large, particularly so when 
the surface tension is negligible. Therefore HySs are pos- 
sible candidates for young millisecond pulsars. 



II. EQUATION OF STATE 

We construct the Equation of State (EOS) of matter at 
high density modeling the Hadronic phase by a relativis- 
tic non-linear Walecka type model |l4| with the inclusion 
of Hyperons [15J. Concerning quark matter, we use an 
MIT-bag like EOS in which quarks can pair to form a 
Color Flavor Locked (CFL) phase El Eft [l^. CFL is 
considered to be the energetically favored type of pair- 
ing pattern, at large densities, in the case of /3-stable, 
electric-charge neutral quark matter. At lower densi- 
ties, a two-flavor Color Superconducting (2SC) phase can 
form, with a smaller energy gap [19j. 

Two scenarios are possible, depending on the value of 
the critical density separating hadronic matter from MP. 
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In the first case the critical density is large, therefore hy- 
perons start being produced at a density smaller than 
the critical one. We will call this phase "hyperon-quark" 
MP [23|. This is the first case discussed below and it 
corresponds to the upper panel of Fig. ^ In the second 
scenario, the MP starts at a density lower than the hyper- 
onic threshold and the hyperon density is much smaller 
than in the previous case. We will explore in particular 
the situation in which hyperons are completely eaten-up 
by quarks and they do not appear in the "nucleon-quark" 
MP (see lower panel of Fig. ^) . In this second scenario 
we will assume that in the MP the density is so low that 
the CFL pairing cannot form and, for simplicity, we will 
consider unpaired quarks in the MP, since anyway the 
bulk viscosity of the 2SC quark phase is similar to the 
viscosity of unpaired quark matter |13| . In both scenar- 
ios, we describe the composition of the MP assuming a 
first order transition and imposing Gibbs conditions |30j . 
Beta stability and charge neutrality are satisfied in every 
density region. 



III. BULK VISCOSITY 

Bulk viscosity is the dissipative process in which a per- 
turbation of the pressure of a fluid element is converted 
to heat. A small variation of the pressure of the system 
can be treated as a perturbation on the densities of the 
different species of particles bringing the system out of 
/^-equilibrium. The reactions between the different par- 
ticles drive the system back to an equilibrium configu- 
ration, with a delay which depends on the characteristic 
time scale of the interactions. Following the formalism of 
Ref. H , we classify all the reactions either as fast or slow. 
Slow processes produce bulk viscosity, because their time 
scale is comparable with the period of the perturbation. 
Fast processes, instead, put constraints on the variation 
of the densities of the different particles. In the follow- 
ing we generalize the formalism of Ref. |8| in order to 
compute the bulk viscosity of MP. The equations needed 
in the computation are rather different in the two sce- 
narios discussed above and we will have to deal with all 
the different cases separately. Moreover, if the effect of a 
non- vanishing surface tension is taken into account, the 
formalism needed to compute the bulk viscosity has to 
be further modified. It is well known that a finite value 
for the surface tension allows the formation of finite size 
structures whose geometrical shape is determined by the 
interplay of the various contribution to the energy, in- 
cluding the Coulomb term |2lj, |22| . A precise estimate of 
surface tension a is unfortunately still lacking, and val- 
ues ranging from a few MeV/fm 2 to a few tens MeV/fm 2 
have been discussed in the literature. We can roughly 
divide this range of values in three windows. Values 
larger than ~ 30 MeV/fm 2 would not allow the formation 
of the mixed phase since it would not be energetically 
favored [TtI l22j . Our analysis is therefore restricted to 
smaller values. In the following we will discuss both the 



case in which the surface tension is so small that it can 
be neglected and the case in which it is non- vanishing. 

A. Negligible surface tension 

Let us first discuss the case in which the surface ten- 
sion is so small that the perturbations of the star, associ- 
ated e.g. to r-modes, can break the finite-size structures 
present in the MP. We have estimated that this case cor- 
responds to values of the surface tension a <1 MeV/fm 2 , 
since to r-modes excitations is associated an energy per 
baryon of a few MeV. 

1. First scenario 

We start our analysis from the first scenario discussed 
in Sec. II, namely from the case in which hyperons are 
present in the MP. We will discuss in particular a param- 
eter set in which only A and E - particles are produced 
in the MP. The only slow processes genera ting viscosity 
are non-leptonic reactions between hadrons [3l| , since the 
largegap prevents weak reactions between quarks. As in 
Ref. (3, we consider the reactions 

n + n < i i^+ p + £~ , (1) 
n+p ^ p + K. (2) 

Following || , we have not taken into account the reaction 

n + n n + A, since the corresponding reaction rate 
cannot be easily estimated. In Ref. 7] it has been argued 
that this rate can be one order of magnitude larger that 
the one associated with reaction (J2J. In the first scenario, 
our results have therefore to be considered as upper limits 
for the viscosity, as it will be clarified in the following. 
Notice anyway that, concerning the damping of r-modes 
instabilities, to be discussed in the following, the most 
important region of the star corresponds to low-moderate 
densities, while the A is produced at larger densities as 
shown in Figs. and 

Concerning fast processes, they come from the follow- 
ing reactions mediated by the strong interaction 

n + A & p + S", (3) 

A + A 3> 2{uds) C FL. (4) 

The last process describes the only possible "melting" 
of hadrons into CFL phase, since the pairing forces the 
number of up, down and strange quarks to be equal. This 
process is fast because of the vanishing value of a. As we 
will see later the "melting" process is instead forbidden 
if the value of a is not negligible. Concerning mechani- 
cal equilibrium, elastic scattering due to strong interac- 
tions, as well as melting processes like the one described 
in Eq.Q), are responsible for (rapid) momentum trans- 
fer between the two phases. The mechanical equilibrium 
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FIG. 1: Particle abundances, as function of the to- 
tal baryon density. The upper panel corresponds to 
the case of hyperon-quark MP and the lower panel to 
nucleon-quark MP. 



between the two phases is reached in a time scale t s de- 
termined by strong interaction, t s ~ 10~ 23 s. The full 
system is, instead, out of mechanical equilibrium on a 
time scale of the order of the period of the perturba- 
tion t p ~ 10~ 3 s. Therefore, during a fluctuation the 
two components of the fluid remain in mutual mechani- 
cal equilibrium. The variations of the pressure in the two 
phases have therefore to satisfy the constraint: 



5P h = SPcFL ■ 



(5) 



The variations Spi of the densities of the various par- 
ticles and the variation Sx of the quark fraction are con- 
strained by the following linearized equations 

= (1 - x)(Sp n + Sp P + Spa + SpT,) 

+ X($Pq) + $X(Pq - Pn ~ Pp ~ PA ~ Pt) , (6) 

= (l-x)(tiPp-5p-£)-6x(p P -p-£), (7) 

= ^2p H Sp H -Pq5p q , (8) 
{H} 

= (3 n 5p n + ppSp p + /3 A Sp A + P^SpT, , (9) 
= aAnSp n + a Ap 5p p + a AA (5pA 

+ CHAT,5pT, - 3aq q 5pq , (10) 



FIG. 2: Hyperon superfluid gap at zero temperature 
and for two values of the baryon density, (see Ref. 0) 



where 



Oiij = U)p, >„.>,, , 
Pi = dP/dpi. 



(11) 
(12) 
(13) 



Eqs. (6), (7) impose baryon number conservation and elec- 
tric charge neutrality. Eq. (8) (where the sum runs on all 
hadrons) imposes the mechanical equilibrium defined by 
Eq. (5). Finally, Eqs. (9), (10) describe the equilibrium re- 
spect to the two strong processes of Eqs. (3), (4). Notice 
that Sx does not appear in Eqs. (8)-(10), because nei- 
ther the pressure nor the chemical potentials explicitely 
depends on the quark volume fraction x- Solving the sys- 
tem allows to express all the dpi and Sx as function of 
Sp n M 

The relaxation time r associated to the weak processes 
reads 



1 



dp 



Sp Sp J Sp 



(14) 



Here Ta and Ts are the rates of the weak interactions 
H3. They have been calculated in Eqs. (4.21), (4.28), 
(4.29) of Ref.Q, where it has been shown that the main 
dependence of the rate on the temperature is a quadratic 
one. Moreover, if hyperon superfluid gaps as the ones dis- 
played in Fig. do develop, they exponentially suppress 
the rates as it will be discussed later. 

The chemical potential unbalance dp is given by 



5p = Sp n - Spa = 25p n - 5p p - 5pT. , 



(15) 



where the variations of the chemical potentials can be 
expressed in terms of dpi as Spi = J2j a ij$Pj ■ A unique 
Sp appears in Eq. (15), since the variations Sp\ and Sp^ 
are constrained by the fast process given in Eq. (3). The 
real part of bulk viscosity, which is the relevant quantity 
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FIG. 3: Bulk viscosity in the first scenario as a func- 
tion of baryon density for various temperatures. Solid 
lines refer to MP viscosity, while dashed lines corre- 
spond to pure hyperonic matter (vanishing hyperon 
gap in the upper panel and finite hyperon gap in the 
lower panel). The dotted lines corresponds to the 
case in which a large hadron-quark surface tension (10 
MeV/fm 2 < a < 30 MeV/fm 2 ) is taken into account 
and they are computed for T= 10 9 ' 5 K (see text). 



relaxation time and as consequence the viscosity too. In 
this sense our results for the viscosity in the first scenario 
must be considered as upper limits. 

We can observe from Fig. [3| that, in the "hyperon- 
quark" MP scenario discussed so far, MP bulk viscosity 
is comparable to bulk viscosity of purely hyperonic mat- 
ter even though the density of hyperons is very small in 
the MP (see upper panel of Fig. ^) [2J| . Notice that in 
the small window of baryonic densities in which both E 
and A hyperons are present, the value of bulk viscosity 
of the MP and of the pure hadronic phase are essentially 
identical. On the other hand, in the baryonic density win- 
p B [f m - 3 ] dows in which only one species of hyperons are present 
in the MP, the bulk viscosity of MP is larger than the 
bulk viscosity of the pure hadronic phase because only 
one decay channel is open (E decay of Eq. or A decay 
of Eq. @) and therefore the relaxation time is larger. It 
is interesting to compute the bulk viscosity in this sce- 
nario also including the effect of hyperon superfluidity, 
as done in Ref. 0. If the energy gaps Ah associated 
with the hyperons are non vanishing, the decay rates are 
suppressed by the factor e~ Aff / T where the gaps have a 
typical shape shown in Fig. For low temperatures the 
viscosity displayes the characteristic features shown in 
the lower panel of Fig. [3J while these features disappear 
at larger temperatures or for vanishing Ah- Concern- 
p B [fur 3 ] ing the appreciable difference between the viscosity of 
MP and of pure hadronic matter for T = 10 9 5 K, this is 
due to the Fermi momentum dependence of Ah , which 
implies that the gaps are suppressed for large hyperon 
densities. In the MP, the density of hyperons is lower 
than in pure hadronic matter and therefore the effect of 
the gaps shows out more dramatically. 

Let us briefly discuss another possible source of bulk 
viscosity related to the formation of a boson condensate. 
Since all quarks in the CFL phase are gapped, the low 
energy excitations are the Nambu-Goldstone bosons as- 



for damping r-modes instabilities, finally reads 
P{loo - 7o)t 



c 



1 + {uJTf 



(16) 



sociated with the spontaneous symmetry breaking of the 
global symmetries In particular, a condensate of tt~ 
or K~ can appear, in the MP, as proposed in Ref. [l8j . 
The role played by these boson condensates in the cooling 
of a CS has been studied in Ref. 23] and the condensates 
can play a role also in the calculation of the bulk viscos- 
ity of the MP. The weak processes involving these bosons 
are: 



where and 70 are the "infinite" frequency adiabatic 
index and the "zero" frequency adiabatic index and u> is 
the angular velocity of the perturbation. 

It is interesting to remark that there are two asymp- 
totic behaviors of the viscosity as a function of the fre- 
quency lo. In the high frequency case {lot 3> 1), the 
viscosity scales as 1/r, while in the low frequency limit 
the viscosity is proportional to r. In the parameter and 
temperature ranges explored in our paper it results that 
we always remain in the low frequency limit, as it can 
be seen from Fig. |3 noticing that the viscosity decreases 
with the temperature. In this regime the addition of an- 
other weak decay (e.g. n + n n + A) decreases the 



K~ 



(17) 
(18) 



In the absence of these bosons, the variation of the den- 
sity of electrons Sp e can be neglected at low temperature 
(T < 10 10 K) because, as already remarked, all leptonic 
reaction rates are much smaller than those associated 
with non- leptonic reactions. If, on the other hand, the 
decay channels of Eqs. (|17fl . (|18|) are open, then the cor- 
responding decay rates are of the same order of the rates 
of the non-leptonic processes of Eqs. Q,© and, there- 
fore, a new contribution to the viscosity appears. While 
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FIG. 4: Bulk viscosity as a function of baryon den- 
sity, for various temperatures in the second scenario. 
Solid lines refer to MP viscosity, while dashed lines 
correspond to pure quark matter. The dotted line cor- 
responds to the case in which a large nucleon-quark 
surface tension (10 MeV/fm 2 < a < 30 MeV/fm 2 ) is 
taken into account and it is computed for T= 10 9,5 K 
(see text). 

a detailed calculation is clearly needed, the main result 
of our work, namely the existence of a large viscosity for 
T < 10 10 K, would even be strengthened. 



P [ frrT 3 ] 




2. Second scenario 

We now discuss the second scenario in which a 
"nucleon-quark" MP forms, hyperons are absent and 
CFL quark pairing cannot take place in the MP [3^ . We 
will assume that in the pure quark matter phase CFL 
gaps can form so that this phase will not contribute to 
the viscosity. The only source of viscosity, neglecting as 
before semi-leptonic reactions, is therefore 



d + u 



Hw 



(19) 



Concerning fast processes, they are given by the "melt- 
ing" of nucleons into unpaired quarks. The chemical 
equilibrium respect to these reactions must therefore be 
satisfied during the perturbation. The linearized equa- 
tions governing the density fluctuations, analogous to 
Eqs. (6)-(10) of the first scenario, read 

= (l-x)($Pn + o~Pp)+x(Spu + dpd + $Ps)/3 

+ $X ((Pu + Pd+ Ps)/3 - Pn - Pp) , (20) 

= (1 - x)$Pp +X(25pu ~ Spa - Sp s )/3 

+ 5x((2 Pu - Pd - p s )/3- pp) , (21) 

= ^2p H Sp H - ^2pqSp q , (22) 

{H} {Q} 

= a pn 6p n + a pp Sp p - 2a uu 5p u - a d dSpd , (23) 

= a nn 5p n + a np 8p p - a uu Sp u - 2a d dSpd ■ (24) 



FIG. 5: Particles densities profiles inside hybrid stars 
of mass M — 1.46M© in the two discussed scenarios. 
Parameters as in Fig. 1. 



Eqs. (20)-(22) impose baryon number conservation, 
electric charge neutrality and mechanical equilibrium as 
in the first case. Eqs. (23), (24) describe the chemical 
equilibrium respect to the two fast processes of "melt- 
ing" . As before we can calculate the chemical imbalance 
associated to the weak reaction (19) and the correspond- 
ing relaxation time. The rate of the process (19) is taken 
from Eqs. (5)-(7) of Ref . [13 (see also footnote j3fj|)- The 
resulting bulk viscosity is shown in Fig. [!]] (lower panel). 
The bulk viscosity of "nucleon-quark" MP is of the same 
order of magnitude of bulk viscosity of strange matter 



B. Effects of the surface tension 



Let us shortly discuss the effect of a non-vanishing 
surface tension a at the interface between hadronic and 
quark matter. We will assume that a < 30 MeV/fm 2 , 
so that MP can form because it is energetically favored. 
Finite-size structures (drops, rods, slabs) form at differ- 
ent densities, to minimize the energy of the MP. In this 
case the perturbation of the pressure due to r-modes is 
too weak to induce the "melting" process which is now 
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a very slow process (in comparison to the period of the 
perturbation) and it plays no role in the calculation of 
the bulk viscosity. The response of the finite size struc- 
tures inside MP to a perturbation of the density corre- 
sponds to three possible processes: a) the formation of 
a drop of "new" phase; b) the merging of two structures 
into a single larger structure; c) the absorption of "old" 
phase into a structure of "new" phase. Obviously, also 
the reverse processes are possible. It is easy to see that 
large values of a suppress all these processes. In partic- 
ular, in case a) the radius of a critical drop of new phase 
increases with a, making it more difficult to produce a 
new drop; case b) shares some similarity with the fis- 
sion problem in nuclear physics, where the process of the 
separation of a heavy nucleus into two lighter nuclei is 
suppressed for larger values of a since, during the fission 
(or merging) process, configurations having a large sur- 
face are produced. Finally, c) can be viewed as a special 
case of b), in which the absorbed hadron can be assim- 
ilated to a small drop of quark matter. Concerning the 
Coulomb interaction, it mainly plays a role in determin- 
ing the size of the structures while it is not so important 
for the response of the structures to the perturbation, 
at least for relatively large values of the surface tension 
a > 10 MeV/fm 2 . In that case, in fact, the screening due 
to electrons almost completely cancels the effect of the 
Coulomb interaction in the nucleation rate |24|. In the 
following for simplicity we will assume that 10 MeV/fm 2 
< a < 30 MeV/fm 2 . In conclusion we reasonably assume 
that large values of a suppress all these processes, there- 
fore the equations describing the melting process as a 
fast reaction, i.e. Eq. (10) (first scenario) and Eqs. (23), 
(24) (second scenario) are not imposed. We can compute 
the viscosity by requiring that the baryon number and 
the electric charge of the two phases are separately con- 
served. This implies that Eq. (6) in the first scenario and 
Eqs. (20), (21) in the second scenario, each of them sep- 
arate into two distinct equations. In the first scenario, 
Eq. (7) is not modified due to the charge neutrality of 
the CFL quark phase. Finally, the equation of mechani- 
cal equilibrium (Eq. (8) in the first scenario and Eq. (22) 
in the second scenario) is still valid and it represents the 
only fast process connecting the two phases. Notice that 
since melting processes are suppressed, the only reaction 
which allows the system to rapidly re-equilibrate is elas- 
tic scattering between the two phases. In the following 
we have assumed that a residual interaction between the 
two phases always exists, similarly to the existence of 
"entrainment" in superfluid neutron matter as discussed 
e.g. in Ref.[2f|. The systems of equations read therefore: 



= 


(1 - x)(5p n - 


f Sp p + 6 p A + 8p^) 






$X(Pn + Pp -+ 


PA + ps) , 


(25) 


= 


\Spq + 5\Pq 


1 


(26) 


= 


(1 - x)(Sp P - 


- Sp s ) - Sx(p P - pt) , 


(27) 


= 


^2p H Sp H - 


PqSpq , 


(28) 



{H} 
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FIG. 6: The chemical unbalance for the second sce- 
nario is shown as a function of the baryonic density. 
The solid line corresponds to the case of a vanishing 
of the surface tension and the dotted line corresponds 
to a finite surface tension. 

= PnSpn+PpSpp + PASpA + fcSpx (29) 
for the first scenario and 

= (l-x)(5pn + 5p p )-5 X (Pn+Pp), (30) 

= x($Pu + $Pd + 8p s ) + 6x(p u + Pd + Ps) , (31) 

= (l-x)Sp P -5 X p P , (32) 

= x(2Sp u - Sp d - Sp s ) + <Sx(2p„ - Pd - Ps) (33) 

= ^PhSph - ^PqSpq (34) 
{H} {Q} 

for the second scenario. 

In Fig. |3 and 01 we show the effect of a non- vanishing 
surface tension on the viscosity (dotted lines). On rather 
general grounds, one can expect that the effect of a non- 
vanishing surface tension is to reduce the viscosity. In- 
deed, the surface tension suppresses the fast processes of 
"melting" which are responsible for the reduction of the 
chemical unbalance. In particular, the only equation con- 
necting the two phases is now the equation corresponding 
to the mechanical equilibrium. Moreover, in the first sce- 
nario the number of constraints on 5pu reduces from four 
to three and, in the second scenario, the constraints on 
SpQ from four to two. It is also interesting to notice that 
in both scenarios, near the beginning of the MP (x — > 0), 
baryon number conservation requires <5x — ^ if the sur- 
face tension is non-vanishing (see Eqs. (26) and (31)). 
On the other hand, in the first scenario pa — at the 
beginning of the MP below the A production threshold, 
and therefore, the constraint 8x — is also satisfied in 
the absence of surface tension, as explained in footnote 
[32]. Therefore the dotted lines coincide with the solid 
lines in both panels of Fig.[3]for ps < 0.6 fm~ 3 . 

In the second scenario, the effect of the surface ten- 
sion is more dramatic. In Fig. we display the chem- 
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ical unbalances corresponding to a vanishing and to a 
finite value of the surface tension, respectively. As al- 
ready remarked Sp/Sp n is larger in presence of a surface 
tension. Concerning the singular behavior of the chemi- 
cal unbalance near the first critical density, it stems from 
neglecting 5p e in the equation of the electric charge con- 
servation. While in general this is a safe approximation 
due to the slowness of the modified Urea process (see 
Ref. in this particular case this approximation im- 
plies the vanishing of the viscosity at threshold. Actually, 
the existence of a finite value for 8p e implies that the vis- 
cosity is small but finite at threshold. It is not possible 
to apply directly the formalism of Ref. Q if two indepen- 
dent perturbations (Sp n and 5p e ) exist in the system. We 
are therefore forced to discuss separately the viscosities 
stemming from the two independent unbalances. Con- 
cerning the viscosity associated with the modified Urea 
process, we remind that it scales as T 6 and therefore it 
is essentially negligible below T ~ I0 10 K while it will 
be included in the calculation of the stability of the star 
presented in the next section. In conclusion, the result 
corresponding to the dotted line of Fig. 0] would not be 
significantly modified by taking into account a finite value 
of Sp e . 



IV. STABILITY OF HYBRID STARS 
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We can now address the problem of the stability of a 
rotating compact star. To compute the critical angular 
velocity we use the standard formalism of Refs.pl 1st l2fij|. 
We need first to integrate the viscosity on the struc- 
ture of the star, which is obtained by solving Tolman- 
Oppenheimer-Volkov equation. In Fig. [5] we show the 
structure of a M = 1.46M Q star for the two scenarios dis- 
cussed above. For simplicity we computed the structure 
of the star assuming a negligible value for a. The main 
effect of the presence of finite size structures is to reduce 
the volume occupied by the MP. For a < 30 MeV/fm 2 , 
which is the limit of validity of our approach, the shrink- 
ing of the MP is rather modest [2l| . The critical angular 
velocity £l C rit is the one for which the imaginary part 
of the r-mode frequency vanishes and it is obtained by 
solving the equation 



1/tgr + 1/tb + 1/t B ([W) = . 



(35) 



Here tgr is the time scale for gravitational waves emis- 
sion while tb and TB(urca) & re the time scales of the 
bulk viscosity produced by hadronic processes and by 
the modified Urea process of the nucleons|3j, respec- 
tively. Results for the critical angular velocity are shown 
in Fig. [7| In the upper panel we compare the stability 
of a purely hyperonic star with the stability of a HyS 
containing hyperon-quark MP. As it can be seen, due to 
the large viscosity of the MP the HyS is as stable as the 
hyperonic star. Let us stress again that our result indi- 
cate that the viscosity in the MP is almost independent 



FIG. 7: Critical angular velocities. The solid lines re- 
fer to hybrid stars, dashed lines correspond to hyper- 
onic stars (upper panel) and to strange stars (lower 
panel). The dotted line corresponds to the case of 
a large hadron-quark surface tension (10 MeV/fm 2 
< a < 30 MeV/fm 2 , see text for details). 



on the hyperonic content and therefore hybrid hyperon- 
quark stars can be stable as long as a tiny fraction of 
hyperon is present. In the lower panel we compare the 
stability of a star made entirely of non-interacting quarks 
with a HyS containing a nucleon-quark MP. We have as- 
sumed that the viscosity of the pure quark matter phase 
in the HyS vanishes, to simulate a CFL core. As it can be 
seen, the small MP region, located near the edge of the 
star, is sufficient to damp the r-modes. For simplicity, we 
have assumed quarks to be unpaired in the MP, but simi- 
lar results should be obtainable if a 2SC phase is present. 
A feature of r-modes is that they are active mostly in the 
outer regions of the star, and therefore the value of the 
bulk viscosity at a not too large densities is crucial for 
the stability of the star. In particular, in the first sce- 
nario the stability of the star is due to the presence of the 
£ hyperons in the pure hadronic phase and in the MP. 
A hyperons, which are produced at larger densities, play 
a lesser role. In the second scenario, a small window of 
MP, present in the outer region, is sufficient to stabilize 
the star. Let us also remark that the star is stable at 
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large temperatures, due to the modified Urea processes 
active in the crust. These processes does not exist in the 
case of pure quark stars, which are therefore unstable 
at large temperatures. Finally, when a finite value for 
the hadron-quark surface tension is taken into account 



the instability window is larger, but the main conclusion 
concerning the stability of a young hybrid star remains 
valid. 

It is a pleasure to thank J.C. Miller and L. Rezzolla 
for very useful discussions. 
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